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Canonical quantization of an action containing curvature squared term requires introduction of 
an auxiliary variable. Boulware etal |l| prescribed a technique to choose such a variable, by taking 
derivative of the action with respect to the highest derivative of the field variable, present in the 
action. It has been shown that pi this technique can even be applied in situations where introduction 
, — , . of auxiliary variables are not at all required, leading to wrong Wheeler-deWitt equation. It has also 

c 2 ^ ' been pointed out that |2| Boulware etal's |l| prescription should be taken up only after removing 

^3 , all possible total derivative terms from the action. Once this is done only a unique description 

Cn ■ of quantum dynamics would emerge. For curvature squared term this technique yields, for the 

first time, a quantum mechanical probability interpretation of quantum cosmology, and an effective 
potential whose extremization leads to Einstein's equation. We conclude that Einstein-Hilbert action 
should essentially be modified by atleast a curvature squared term to get a quantum mechanical 
formulation of quantum cosmology and hence extend our previous work [H| for such an action along 
with a scalar field. 
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I. INTRODUCTION 



Quantum cosmology is most elusive as the role of time [pj is not unique and one is being unable to define the 
^ Hilbert space. This is due to the fact that 'time' is not an external parameter in general theory of relativity, rather 
^ , it is intrinsically contained in the theory, unlike its role in quantum mechanics or quantum field theory in flat 
O^' space time. In curved space time, different slices correspond to different choices of time leading to inequivalent 
; I \ quantum theories. Likewise, the canonical quantization of gravity is devoid of an unique time variable and hence 
bJ[)i the definition of probability of emergence of a particular Universe out of an ensemble is ambiguous. The canonical 
quantization of Einstein-Hilbert action together with some matter fields yields the Wheeler-deWitt equation which 
does not contain time a priori, although it emerges intrinsically through the scale factor of the Universe. However, if 
the canonical variable is so chosen that one of the true degrees of freedom is disentangled from the kinetic part of the 
J^ . canonical variables, then this kinetic part in the corresponding quantum theory yields a quantum mechanical flavour 
of time. This is possible only if the Einstein-Hilbert action is replaced by curvature squared action or modified by the 
introduction of curvature squared term, in the Robertson- Walker minisuperspace model. In a recent publication H 
we have presented such a choice of the canonical variable in a curvature squared action whose quantization yields a 
Schrodinger like equation. Further, as a consequence, we could identify the nature of the space and time like variables 
in the Robertson- Walker minisuperspace model from the equation of continuity and hence could establish the idea of 
the probability and the current densities in the context of quantum cosmology. 

The relevance of the fourth order gravity theory in cosmology was first revealed by Starobinsky B, though it 
appeared earlier in the context of quantum field theory in curved space time. Starobinsky IJ] presented a solution of 
the inflationary scenerio without invoking phase transition in the very early Universe, from a field equation containing 
only the geometric terms. However, the field equations could not be obtained from the action principle, as the terms 
in the field equations are generated from the perturbative quantum field theory. Later, Starobinsky and Schmidt 
Q have shown that the inflationary phase is an attractor in the neighbourhood of the fourth order gravity theory. 
Such wonderful relevance of higher order gravity theory in the context of cosmology inspired some authors to give 
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interpretation of the quantum cosmological wavefunction with R^ |g] and even R^ terms in the Einstein- Hilbert 
action. 

Further the functional integral for the wavefunction of the Universe proposed by Hartle-Hawking Iq] runs into 
serious problem, since the wavefunction diverges badly. There are some prescriptions ||9[ to avoid such divergences, 
though a completely satisfactory result has not yet been obtained. However, to get a convergent functional integral, 
Horowitz |lC|] proposed an action in the form 

S^ j d^X^[AClki+B{R~A\f] (1) 

where Cijki is the Weyl tensor, R is the Ricci scalar, A is the cosmological constant and A, B are the coupling constants. 
The action (1) reduces to the Einstein-Hilbert action at the weak energy limit flG]. To obtain a workable and simplified 
form of the field equations, one may consider a spatially homogeneous and isotropic minisuperspace background, for 
which the Weyl tensor trivially vanishes. In a previous paper B, we considered the action (1) retaining only the 
curvature squared term. The field equations for such an action can be obtained by the standard variational principle. 
In the variational principle, the total derivative terms in the action are extracted and one gets a surface integral which 
is assumed to vanish at the boundary or the action is chosen in such a way that those surface integral terms have 
no contribution. However, for canonical quantization this principle is not of much help and one has to express the 
action in the canonical form, which is achieved only through the introduction of the auxiliary variable. Auxiliary 
variable can be chosen in an adhoc manner, and different choice of such variable would lead to different description of 
quantum dynamics, keeping the classical field equations unchanged. In view of this Ostrogradski pl| in one hand and 
Boulware etal [|l| on the other, have made definite prescriptions to choose such variables. Ostrogradski's prescription 
PI was followed by Schmidt [12|. We shall consider Boulware etal's [|l| prescription in which it was proposed that the 
auxiliary variable should be chosen by taking the first derivative of the action with respect to the highest derivative 
of the field variable present in the action. 

Hawking and Luttrell ||l^ utilized Boulware etal's ||] technique to identify the new variable and showed that the 
Einstein-Hilbert action along with a curvature squared term reduces to the Einstein-Hilbert action coupled to a massive 
scalar field, assuming the conformal factor. Horowitz pOf on the other hand, showed that the canonical quantization 
of the curvature squared action yields an equation which is similar to the Schrodinger equation. Pollock [|l4| also 
used the same technique to the induced theory of gravity and obtained the same type of result as that obtained by 
Horowitz ||10| , in the sense that the corresponding Wheeler-deWitt equation looks similar to the Schrodinger equation. 

The striking feature of Boulware etal's ||l[ prescription is that it can even be applied in situations where introduction 
of the auxiliary variable is not at all required, eg. in the induced theory of gravity, vacuum Einstein-Hilbert action etc. 
It is not difficult to see that the classical field equations remain unchanged with or without the introduction of the 
auxiliary variables. Now the question is whether quantum dynamics remain unaffected? Putting other way round, the 
question turns out to be whether the solution of the Schrodinger like equation obtained by Pollock Q in the induced 
theory of gravity by introducing auxiliary variable will satisfy the Wheeler-deWitt equation corresponding to the same 
action that can be obtained without introducing such variable? The answer is simply no. We have shown that 
the introduction of an auxiliary variable does not change the classical field equations in situation mentioned above, 
however, the quantum dynamics are different. The vacuum Einstein-Hilbert action in the background of Robertson- 
Walker metric has been chosen for this purpose. We have shown in such a toy model H], 'toy', since it is not required 
to follow Boulware etal's |l|| prescription in the vacuum Einstein-Hilbert action, that the introduction of an auxiliary 
variable leads to wrong Wheeler-deWitt equation. So what went wrong with Boulware etal's M prescription? This 
is what we have concluded in |g], that before taking up the said prescription, one has to eliminate all possible total 
derivative terms from the action. Once this is done, it would not at all be possible to introduce auxiliary variables in 
situations mentioned above, where such introduction is unnecessary, and hence only a unique and correct quantum 
dynamics would emerge. In this sense, the auxiliary variable identified by Hawking and Luttrell p3| and Schrodinger 
like equation obtained by Horowitz [Q and Pollock ||lj] are wrong. We g] have obtained the correct quantum 
dynamics for curvature squared action and showed that the continuity equation identifies the nature of space and 
time like variables in the minisuperspace, establishing, perhaps for the first time the quantum mechanical idea of the 
probability and the current densities in quantum cosmology. It has also been shown that the correct description of the 
curvature squared action leads to an effective potential, whose extremum yields vacuum Einstein's equation, which is 
a desirable feature of higher derivative gravity theory, in the weak energy limit. The emergence of the Schrodinger like 
equation, the quantum mechanical concept of the probability density and the effective potential lead us to conclude 
that the quantum description of cosmology is incomplete without atleast a curvature squared term in the action. This 
is the reason why in this paper we have taken up the task of generalizing our previous work O] by modifying the 
Einstein-Hilbert action with a curvature squared term along with a matter field. 

As already mentioned, in order to obtain the correct and unique quantum cosmological description of a system whose 
action contains curvature squared term, one has to first remove all possible total derivative terms from the action and 



then follow Boulware etal's |l| prescription. The whole formulation requires four steps, viz. 1) remove all possible 
total derivative terms from the action, 2) introduce auxiliary variable, which is the first derivative of the action with 
respect to the highest derivative of the field variable present in the action, 3) cast the action in the canonical form 
and obtain classical field equations, 4) introduce basic variables to canonically quantize the Hamiltonian constraint 
equation. At this stage we would like to mention that the Hamiltonian in our formalism before introducing basic 
variables [^ turned out to be the same as that obtained by Schmidt Jf^ following Ostrogradski's ||f^ prescription. 



This is due to the fact that, in that prescription |fl| total derivative terms do not appear. However, the importance 



of removal of the total derivative terms has not been stressed, the method of introducing auxiliary variables and 



formalism and 
Therefore we 



process of casting the action in the canonical form is not so vivid and elegant as in Boulware etal's Jl 
quantization procedure by introducing basic variables has not been mentioned in that approach [[11 
have selected the other method. 

This paper has been organized in the following manner. In section 2. we present the field equation for an action 
containing curvature and curvature squared term along with a minimally coupled matter field, following the proposal 
mentioned above. Further we cast the action in canonical form with respect to the basic variables, which shows the 
appearence of a Hamiltonian —Pa of a system described by the basic variables and their conjugate momenta, though 
the actual Hamiltonian H ie. the super Hamiltonian constraint to vanish. The form of the above action also identifies 
the time variable, which is the intrinsic time, whose canonical momentum is Pq, or the Hamiltonian —Pa can be 
thought of as the generator of the translation along the above time variable. The action written in this form indicates 
the true dynamical degree of freedom, disentangled from the kinetic variable ie. the time variable. We claim that the 
Hamiltonian {—Pa) is the correct Hamiltonian in the sense that it contains the basic variables and their conjugate 
momenta and its canonical quantization yields the Wheeler-deWitt equation which mimics the Schrodinger equation, 
along with an effective potential. 

In section 3. we present the conventional continuity equation of quantum mechanics from the Schrodinger like 
Wheeler-deWitt equation. The continuity equation in quantum cosmology was first presented by Sanyal and Modak 
only in a recent publication and it is interesting to note that, as the probability and current densities are finite, so 
also the wavefunction and its derivatives are finite throughout the evolution of the Universe if and only if there are 
no singularities in the domain of quantum cosmology. 

Section 4. has been devoted to extract physical implication of the effective potential Ve- There are singularities 
in the effective potential. It is the function of the scalar field and the expansion parameter, which is the Hubble 
parameter of the classical field equations and it is diverging both at the zero's and at the very large values of the 
expansion parameter. The effective potential is asymetric with respect to the expansion and collapse of the Universe 
and its extremum yields classical equations. 

Section 5. is devoted in making some concluding remarks. 

II. CLASSICAL FIELD EQUATIONS AND QUANTIZATION 

We consider the following action 

S = j d^x./—g[-^{R - ^R') + ^(^0,^ </>'" - Vim- (2) 

Under the following form of the closed Robertson- Walker metric 

ds^ = e^^^") [dr]^ - dx^ - sin\{de'' + sin^Od^'^)] (3) 

the action (2) can be expressed as 

S^J rf^[^{(l + «' + ")e'" - /3(1 + a^+ar} + ^02^20 _ v{cf,)e'% (4) 

where, d = 4^. According to our proposal, we remove all removable total derivative terms from the action, so that 
only lowest order terms appear in the action (4). Thus, 

S^ f dry[Af {(1 - d2)e2" - /?(1 + a^f - pa^} + i(/.2g2a _ v{<f)e^"] + 5™, (5) 

where, M = |g and s™ = M[ae^°' - 2/3(d + ^)]. 

At this stage, we define the auxiliary variable, which is 



dS 
Mf3Q = - ^— = 2M/3a, ie.Q = 2a. (6) 

oa 

Introducing the auxiliary variable Q in aetion (5) and expressing the action in the canonical form, after removing the 
remaining total derivative terms, we obtain, 

5 = /" dr/[M{/3Qd + (1 - d2)e2" - ^{1 + a^f + '^Q^} + i^^e^" - V{4>)e^''] + S^, , (7) 

where, Smi is the surface integral, given by Sm-^ ~ Sm — PQa. Assuming Smi vanishes at the boundary, the classical 
field equations are found as, 

f3{Q - I2aa^ - 4d) - 2e2"(l + a^ + a) - -jri^'^ - W{cj))e^'^) == 0, (8) 

Q = 2d, (9) 

+ 2d(/) + e2"^=O, (10) 

a(p 

/3{3a^ + 2a^-aQ + ^-l) + (a^ + l)e^" = ^(^^^ ^ v{cf>)e''^), (11) 

the last one being the Haniiltonian constraint equation. The definition of Q given in the equation (6) is thus recovered 
in equation (9). While equation (10) is the well known continuity equation of the scalar field 0, equations (8) and 
(11) are the correct classical field equations as it is evident by replacing Q by 2a. Now the Hamiltonian constraint 
equation in the phase space variables can be expressed as, 

where, Pa, Pq and P^ are the canonical momenta corresponding to a, Q and (j) respectively, and d, Q and cf) are given 

by 

The Wheeler-deWitt equation is obtained through the quantization of the Hamiltonian constraint equation (12). Now 
as already mentioned, canonical quantization should be performed with the basic variables, viz. a and d. Let us 
choose d = x, and hence replace Pq by M(3x, in view of the first equation given in (13). Further, from equation (6), 
we get 

,,^^ dS dL dL 

MpQ^-— = -— = -— = -P,, 14 

oa oa ox 

which means that Q should be replaced by —jf^- From equation (14) it is evident that Px is the canonical momentum 
corresponding to x. Hence equation (12) turns out to be 

^ = a^^a - ^ + ^P|e-2" + M/3(l + x'f + M{x' - l)e2" + y(0)e4". (15) 

As the Hamiltonian H is constraint to vanish, we get 

AMpx 2x ^ X X X 

We have already mentioned that P^, P^ and Px are the canonical momenta corresponding to a, 4> and x, respectively, 
however to understand the physical meaning of Pq, we write the action (7) in terms of the canonical coordinates, 
apart from a surface term, remembering iJ = 0, as 



S= f Ldrj= j[P^a + PqQ + P^</.]d?7 (17) 

Further, from equation (14) we find 



1 dP^ a dP^ 



MP dr] MP da 



(18) 



Now using the expressions for Pq and Q given by equations (13) and (18) respectively, the relations da — adrj, (p — -£^a 
in action (17) and finally rearranging terms, we obtain 

Apart from the surface term the action (19) looks similar to the conventional form of the action of a mechanical 
system described by the configuration space variables x, ^ and their canonical conjugate momenta Px,P^, with a 
Hamiltonian (—Pa), which is also the canonical momentum of a parametrized by the variable a, though the true 
Hamiltonian H of the gravitational system constraint to vanish. We know that the Hamiltonian of a system is the 
generator of the time translation. Since, — Pq in our system plays the role of the Hamiltonian, hence it acts as the 
generator of translation along a, and the variable a acts as the time variable. This feature is also consistent with 
the intrinsic concept of general theory of relativity, as time has no independent existence from geometry in describing 
gravitation, rather it is inbuilt in the theory, unlike situations encountered in the conventional classical and quantum 
mechanics, where time is an external parameter. The functional form of the Hamiltonian {—Pa) is given by (16) and 
it is a function of x, cj), Px, P^ and the intrinsic time parameter a. The canonical quanization of (16) yields 

where n is the operator ordering index and 

K = K(x,0,a) = ^(.^ + l)^ + ^(.^-l)e- + ^Me-. (21) 

XX X 

Equation (20) is the correct Wheeler-deWitt equation corresponding to an action, containinng curvature and curvature 
squared terms along with a minimally coupled scalar field. It looks similar to the Schrodinger equation, though it 
differs from that obtained by Pollock |14| . 

III. PROBABILITY AND CURRENT DENSITY 

One of the most important feature observed in quantum mechanics is that, the state of a system is described by a 
wavevector belonging to an abstract Hilbert space and the norm of the wavevector must be positive definite or zero. 
This idea emerged from the interpretation of the probability density to describe a given state from the continuity 
equation which is obtained by using the Schrodinger equation. Probability interpretation follows from the simple 
mathematical appearence of the Schrodinger equation. No such interpretation of the probability density in general 
exists in quantum cosmology, when the action contains terms linear in the Ricci scalar coupled with some matter 
field. This is due to the fact that there is no time a priori in the Wheeler-deWitt equation in a gravitational theory 
described by the action (2) with /3 = 0. It is to be noted that the continuity equation along with the conventional 
notion of the probability density can only be introduced with the proper choice of the canonical variables in such a 
way that the Hamiltonian constraint is quadratic in the canonical momenta along with a term linear in momentum 
whose canonical coordinate acts as an intrinsic time variable. This type of canonical quantization is possible only 
when the action contains atleast a quadratic curvature term. 

Equation (20) can be writtn as 

ih^ = Ho^J, (22) 



where 



h^ , &■" n d . h^ _2a d'^ 



""^^mrx^dx^^xd-J-Tx^ w^'^'^'^'^^''^- ^''^ 



and 



It is to be noted that Hq operator is hermitian and as a consequence we obtain the continuity equation, viz. 

|^+V.J = 0, (24) 

where p and J are the probabihty and the current densities respectively for the choice of the operator ordering index 
n = —1. p and J are given by p = ip*ip and J = {Jx, J^, 0), where 

'^- = ^(^*^--V'V':), (25) 

^^^-^^c^v^-v-^;), (26) 

One can also find the continuity equation for other values of the operator ordering index but with respect to a 
new variable y which is functionally related to x only. It should be mentioned here that the continuity equation in 
quantum cosmology was presented by Sanyal and Modak for the first time in a recent publication. In the present 
work we have obtained the continuity equation once again. Since the above probability and the current densities are 
dynamical quantities, therefore the wavefunction and its derivatives should remain finite at all epoch of the evolution 
of the Universe, if and only if there are no singularities in the domain of quantum cosmology. 

Further in analogy with the quantum mechanics it is noted that the variable a in equation (20) can be identified 
as the time parameter in quantum cosmology. The variables x{— a) and (f) act as spatial coordinate variables in 
the context of quantum cosmology. The coordinate variable x{— a), which is just the expansion parameter and the 
Hamiltonian operator Hq diverges at the bounce of the Univese, ie. at d = 0. 

IV. EFFECTIVE POTENTIAL AND ITS EXTREMIZATION 

The effective potential given by equation (21) is found to be assymetric with respect to the expansion parameter 
(x), ie. Ve{—x) — —Vei+x) and it is also unstable both at the short range and the long range values of x.The 
significance of the effective potential becomes clear when it is extremized, yielding interesting results. In the weak 
energy limit, ie. when the dcBroglie wavelength is much larger in comparison with the average separation (ie. far 
below the Planck scale) the contribution of the kinetic energy is sufficiently small with respect to the potential energy 
in the Hamiltonian . At this regime the Hamiltonian is almost dominated by the potential energy. Now the extremum 
of the effective potential Ve given by equation (21), with respect to x yields, keeping cj) and the time parameter a 
fixed, 

(3{3x^ + 2x^ - 1) + (a;2 + l)e^°' - ^V{<j)) = 0. (28) 

If V{(j)) has got no extremum then in the limit f3 —>■ 0, ie. in the absence of the curvature squared term, equation (28) 
takes the following form 

x^ + l^a^ + l = ^V{^), (29) 

which is exactly the classical constraint equation in the Einstein's gravity, along with a scalar field, apart from a 
kinetic energy part corresponding to the scalar field, that we have already assumed to be small. Further, for the 
action dominated by the curvature squared term, ie. at /3 ^ oo, we have 

a;^ + 1 = 0, (30) 



x^ 



l-O. (31) 



These pair of equations have already been obtained in 0. Here, we have just been able to recover them in the limit 
P — > oo, ie. at the epoch of the evolution of the Universe dominated by the curvature squared term. These equations 
imply that the extremum of the effective potential Ve is at a coordinate position x — a, satisfying either equation 
(30), which is the vacuum Einstein's equation that admits Euclidean wormhole solution, or equation (31), which is 
the epoch at which the Universe admits a solution, given by a = ^ ~^' , 'a' being the scale factor in proper time 't'. 
For this solution the horizon radious r^j is proportionaql to ln(t), and so as t ^ 0,r/f — > cxd, and thus the horizon 
problem is solved. The effective potential Ve diverges in both the directions of x — > and x — > cx), even for /3 — > cx). 
Further, it can be shown that ^^^ > 0, for condition (31), ie. Ve has got a minimum. These altogether imply that the 
Universe, while sitting at the minimum of the effective potential, will expand at the rate a proportional to t, solving 
the horizon problem. 

Now, had there been an extremum of the scalar field potential V{4>), then extrcmizing the effective potential Ve, 
given in equation (21), with respect to </>, keeping x, a fixed, we get 

e^" dV , , 

Mx del) ^ ' 

The condition that a function (Ve) of two variables x, and (j> (not considering the time parameter a) has got a 
minimum, is 

d^Ve d'^Ve d'^Ve 

dx"^ dcj)^ dxdcj) 
The left hand side in the present context turns out to be 

|^[e-+2/3(3.g + l)]^U.,.,„ (34) 

which is positive definite for ^ A \at4>=<po > 0^ i^- ^r a minimum of the scalar field potential. Therefore the effective 
potential has got an extremum located at the configuration space variables (fi — (Jjq and x = xq, determined by equation 
(28). 

It can further be shown that g^f \<p.a=constant > 0, which means that the effective potential is truely a minimum at 
the location of the configuration space variables (f) — (J)q, x — xq, determined by equation (28), where V{(j)) — V{(J)q). 
Further, if we choose the minimum of the scalar field potential V{(t>Q) — 0, then equation (28) reduces to either 

a;2 + 1 = (35) 



/3(3a;2 - 1) + e^" = 0. (36) 

these pair of equations again imply that the Universe is sitting at the minimum of the effective potential 14, at the 
locations of the configuration space variables 0o and x = xq, satisfying either equation (35), which is the vacuum 
Einstein's equation admitting Euclidean wormhole solution, or equation (36), that admits an oscillatory solution, in 
the form 

a^^Ps^n\'-^\ (37) 

where, as already mentioned, 'a' is the scale factor in proper time 't'. 

V. CONCLUDING REMARKS 

In summary, we have observed that, different choice of the auxiliary variable, leads to different inequivalent quantum 
dynamics, keeping the classical field equations unchanged. Hence, we have suggested a principle of choosing the 
correct auxiliary variable and the true degree of freedom, which would lead to the correct and unique quantum 
description. For this we have modified Boulware etal's principle M by adding only one statement at the begining 
that, before choosing such variables one has to remove all removable total derivative terms from the action. That 
this would lead to the correct quantum description, has been proved in a toy model, in a recent publication by 



Sanyal and Modak g . In addition, it has been shown that the quantum dynamics obtained by this principle leads to 
certain desirable features viz. quantum mechanical probability and current density interpretation from the continuity 
equation of quantum cosmology and an effective potential whose extremum yields classical field equations. It seems 
that a quantum mechanical probability interpretation is possible, only if the Einstein-Hilbert action in the Robertson- 
Walker minisuperspace model is modified by atleast curvature squared term, which appears in the one loop correction 
of perturbative quantum field theory in curved space time, and is the most dominant term in the quantum domain. 
Thus we may conclude that, quantum mechanical probability interpretation of quantum cosmology is a generic feature 
of curvature squared gravity, at least in the Robertson- Walker minisuperspace model. 
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